Limits to superweak amplification of beam shifts 
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The magnitudes of beam shifts (Goos-Hanchen and Imbert-Fedorov, spatial and angular) are greatly enhanced 
when a reflected light beam is postsclcctcd by an analyzer, by analogy with superweak measurements in 
quantum theory. Particularly strong enhancements can be expected close to angles at which no light is 
transmitted for a fixed initial and final polarizations. We derive a formula for the angular and spatial shifts 
at such angles (which includes the Brewster angle), and we show that their maximum size is limited by 
higher-order terms from the reflection coefficients occurring in the Artmann shift formula. © 2013 Optical 
Society of America 
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The law of specular reflection fails for optical beams, 
as any physical beam is made up of multiple plane wave 
components, each experiencing a different reflection co- 
efficient. This variation leads to small deviations from 
perfect specular reflection: an angular shift in the beam's 
propagation direction [1], and a spatial shift of the cen- 
ter of the reflected beam. These beam shifts have at- 
tracted much experimental and theoretical interest in re- 
cent years (reviewed in [2]), consisting of shifts in the in- 
cidence plane (Goos-Hanchcn shift), perpendicular to it 
(Imbert-Fedorov shift), and of various modifications for 
different beam profiles (whilst assuming the incident in- 
tensity is axisymmetric) . In appropriate units, the shifts 
are independent of the particular radial beam profile, 
provided the beam is paraxial. 

As a wave optical phenomenon, beam shifts arc typi- 
cally small (a few wavelengths for spatial shifts, and ap- 
proximately the beam's Fourier width for angular shifts) . 
They are largest for incidence close to the critical an- 
gle or the Brewster angle, due to the special form of the 
Fresnel coefficients at these angles [3,4]. Here we describe 
another way of enhancing the magnitude of beam shifts, 
namely by 'postselecting' a particular polarization state 
with an appropriate choice of the polarization analyzer. 
The payoff for the large apparent shift in the chosen po- 
larization output is a small integrated intensity in that 
component; this consequence is standard in the quan- 
tum mechanical theory of 'weak measurements' [5—7] , to 
which beam shifts are a classical analog [8,9]. 

In this spirit, we define a null-reflection angle with re- 
spect to a chosen initial and final polarization state, to be 
an angle of incidence for a reflected plane wave to be zero 
after the analyzer. Such angles are somewhat analogous 
to the Brewster angle at which plane waves polarized 
in the plane of incidence (p-polarized) are fully transmit- 
ted: the corresponding reflection coefficient rp(0B) = 0, 
giving zero refiected light projected onto an analyzer set 
in the p-polarization state [10]. Null-reflection angles are 
not polarizing angles [11, 12], but they are interesting 
in the context of beam shifts, as they occur for arbi- 



trary states of elliptic polarization, when the analyzer is 
crossed with (orthogonal to) the reflected polarization, 
and they can be found in regimes of both total and par- 
tial reflection. 

Beams, made up of many plane wave components, do 
not have a homogenous polarization after reflection, and 
hence the projection onto an analyzer polarization is 
nonzero but possibly shifted. With appropriate choice 
of initial and final polarizations, the beam shift close to 
Brewster or null-refiection angles can be very large. In 
the quantum analogy, this is the 'superweak' regime [13], 
where the value of the beam shift in the postselected 
component takes on a larger value than the shift of to- 
tal intensity for either s or p incident polarization, albeit 
with a greatly reduced intensity. For general incident and 
analyzer polarizations, the net component shift is arbi- 
trary (it has a spatial and angular shift both with lon- 
gitudinal and transverse components) , and the apparent 
shift of the total beam is an appropriately weighted sum 
of orthogonal analyzer polarizations [8]. 

The magnitudes of optical beam shifts are usually cal- 
culated using the approximate Artmann formula [14] 
suitably generalized to angular, transverse and post- 
selected shifts [8]. At Brewster and null-reflection angles, 
this reflection coefiicicnt is zero, meaning the Artmann 
formula blows up, implying that the shift is arbitrarily 
large. A more refined calculation, given below, includes 
the next order in the Taylor expansion of the reflection 
coeflicient, regularizing the singularity to be a resonance- 
like curve [15, 16] whose maxima and minima give the 
largest values of the beam shift, which now depend on 
the shape of the beam. This behaviour is not unique to 
the Brewster angle, and our results extend previous the- 
oretical work on Brewster refiection [17, 18]. 

For plane waves initial polarization and post-selecting 
analyzer are represented by Jones vectors E and F re- 
spectively, and the angle of incidence is 6. The reflection 
is described by a 2 x 2 Jones matrix R, which, in the 

p,s basis is given by R = o" ^s) > where rp — rp{9) 

Vg = rs{6) are the Fresnel reflection coefiicients [10]. A 
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reflected plane wave will therefore have a polarization 
R- E, and the amplitude after the analyzer is i^* • R • 
At ^B, this is zero since E is Bp and J'pC^b) = 0, giving 
zero transmission regardless of F. At null-reflection an- 
gles, F is orthogonal to the reflected polarization R • E, 
and so any incident can be made into such an angle 
with an appropriate choice of F given E. 

A reflected beam is represented by a bundle of plane 
waves, labelled by transverse wavevectors K, centred on 
a central plane wave with wavevector K ~ in our cho- 
sen beam coordinate system. Following the approach of 
Ref. [8] , it is convenient to work in cylindrical coordinates 
with axis the 'virtual reflected beam' (i.e. unshifted), 
with azimuth (p and radius r in real space, azimuth a 
and radius sin 6 in Fourier space. The shift formulas for 
a reflected beam follow from a Taylor expansion of 



p{K) = F* ■ RiK) ■ E, 



(1) 



where R(-K') is the fC-dependent reflection matrix, 
which is a 2 X 2 Jones matrix in beam coordinates, 
and R(0) is the reflection matrix of the central plane 
wave. Since the beam is paraxial, its spectrum is local- 
ized around the origin in Fourier space, justifying the 
expansion of p{K) up to only the flrst few terms [2], 



p{K) « po+Spi- 
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, (2) 



where c(a) — cosa,s{a) = sin a, and po,Pi and P2 all 
originate from the Taylor expansion of p. We will implic- 
itly use the ep, Cg basis, with polar coordinates r, cf) in real 
space and K = K{kS, a) in Fourier space (with the range 
of 5 so small that sin J S). The unperturbed, virtual 
amplitude (which corresponds to the incident beam) is 
real and axisymmetric, and is written ip(r), with Fourier 
transform <f{K). We assume that these are normalized, 
^^r\^{rrdr=J^K\^{KrdK^l. 

The beam shift is determined by the centroid (posi- 
tion expectation value) of the reflected beam il^{r), which 
is no longer axisymmetric. In Fourier space, p{K) is a 
multiplication operator on the amplitude, and i^{K) = 
p{K)(p{K). A good approximation of the angular shift 
A, valid for a paraxial beam close to the Brewster and 
null-reflection angle, comes from applying the expansion 
(2) to this expectation, i.e. 



l J^^dKKj^''daK\jiK)\^ 
k j;;^dKKj^'''damK)\^ ^ 
Jf^ dS J^^ da (5^(cosQ;,sinQ;)|-0(lir)'^ 



(3) 
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In this expression we have omitted small terms such 
as those combining po with third order derivatives of 
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Fig. 1. Angular shifts close to a conventional Brew- 
ster angle. The 6'-dependent angular shift (thick, red 
curve), with = (1, 0) is compared with the shift of the 
total intensity [1] (light pink curve) for a reduced refrac- 
tive index of n = 1.5 (e.g. glass-air interface). Both shifts 
follow a resonance-like curve close to the Brewster angle 
at 6b = 56.31°. Black circles indicate obtained by 
numerical integration of (3). The moments of the spec- 
trum ((5^) and ((5"*) are obtained for a Gaussian beam 
with beam waist of 50 pm and a wavelength of 0.526 
pm. On the same abscissa, we plot the ratio of reflected 
to incident intensity (thin, gray curve), showing a mini- 
mum close to the maximum shift. The flgure also shows 
three contour plots (insets) of the angular spectrum of 
the light beam at angles 6 = 56.02° (a), 56.32° (b) and 
57.87° (c). Orange crosses indicate the origin from which 
the shifts are measured. The strength of the contour is 
shown in the colorbar in the top right corner. 



p. The usual Artmann formula for this case contains 
only terms with po and pi. In the last line, (5") = 
fc^ /o°° dS (5""'"^|i^((5)p, the nth radial moment of the spec- 
trum, and ((5"+^) ^ {S") for each n, normalized such 
that {S'-^) = 1. Thus, in usual beam shift circumstances, 
only the flrst term in the numerator and denominator 
contributes to the angular shift, giving the known angu- 
lar shift formula [1], which is independent of the beam 
proflle apart from the overall scaling by the second mo- 
ment of the spectrum. However, when |po| <^ 1, close to 
the Brewster and null-reflection angle, the higher terms 
included in Eq. (4) contribute, replacing the singular- 
ity with a resonance-like complex singularity, as shown 
in Fig. 1. The postselection enhances the already large 
Brewster angle shift for total intensity, which follows a 
similar curve [18]. Beam shifts at null-reflection angles 
can be both angular and spatial, and in an arbitrary di- 
rection in the xy-plane [8]; the shift is purely angular, in 
the plane of incidence at real Brewster angle. 

The calculation of the formula for the spatial shift D 
at the null-reflection angle is similar to the angular calcu- 
lation above, but with ip(r) replacing its Fourier trans- 
form. However, in real space, the reflection operator p 
becomes a differential operator with K — >■ — iV in the 
expansion (2), giving ip{r) = p{—i'V)ip{r). The action of 
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the operator in Taylor expanded form on ip{r) is easy to 
calculate, and so the formula for the regularized spatial 
shift D is 



D 



(5) 



1 Im |poPi 
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P*itrP2)} 
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Re(p5trP2) 



(6) 



omitting the technical but straightforward integrals and 
truncating at second order derivatives of p. Unlike the 
usual beam shift, the regularized spatial shift now de- 
pends on the beam's Fourier width {S'^). As before, when 
p is not small, D is given by Artmann formula for the 
spatial shift [14]. However, when |po| ^ 1, once again the 
next terms in numerator and denominator play a role, 
regularizing the infinite expression in the usual Artmann 
formula. 

To illustrate the amplification in the superweak 
regime, in Fig. 2 we plot the post-selected Imbert- 
Fedorov shift around a null-reflection angle for initial 
circular and final 45° linear polarizations, which again 
follows a resonance-like curve. This shift is much larger 
than the largest possible Imbert-Fedorov shift in the to- 
tal intensity at this angle, which occurs for eigenpolariza- 
tions of the Imbert-Fedorov effect [19]. Because the weak, 
postselected shift is larger than this shift, we truly are in 
the superweak regime, where weak values are far larger 
than the spectrum of the operator. The concurrent low 
intensity, however, is still experimentally accessible [20]. 
We note that this choice of polarizations leads to an op- 
tical vortex in the reflected beam. 

For usual beam shifts, the spatial and angular shifts 
are the real and imaginary part of the same complex ex- 
pression, which may be interpreted as a complex weak 
value of the Artmann operator [9,21,22]. However, this 
is no longer the case for a general beam profile, as the 
weighting on the spectral moments is different in Eqs. (4) 
and (6). This discussion illustrates another fact about 
shifts close to the Brewster angle: in going beyond first 
order, the detailed behavior of the shift becomes depen- 
dent on the spectral profile, unlike the regular beam 
shift. We note that this is not necessarily the case for 
higher-order beam shift phenomena, as the case of a 
high-order optical vortex on reflection demonstrates [23] . 
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Fig. 2. Large, weak spatial shift close to a null-reflection 
angle 6^1- The shift Dy (thick, red curve) is a function 
of 6 for a refractive index of n = 3/10, showing the 
resonance-like curve around ^nr ~ 41.76°. The light pink 
curve is the shift for Imbert-Fedorov eigenpolarization 
which is smaller than 0.15/Ltm across the shown range. On 
the same axis, the ratio of reflected to incident intensity 
is plotted (thin, gray curve). The black circles are the 
shift for a numerical experiment for a Gaussian beam 
with a waist of 50 fim and a wavelength of 0.526 pm. 
Inset are contour plots of the beam's amplitude at angles 
= 41.60°(a), 41.77°(b) and 42.86°(c). Orange crosses 
indicate the origin for the contour plots. 
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